Symmetry breaking in the Hubbard model at weak coupling 
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The phase diagram of the Hubbard model is studied at weak coupling in two and three spatial 
dimensions. It is shown that the Neel temperature and the order parameter in d — 3 are smaller than 
the Hartree-Fock predictions by a factor of q — 0.2599. For d = 2 we show that the self-consistent 
(sc) perturbation series bears no relevance to the behavior of the exact solution of the Hubbard 
model in the symmetry-broken phase. We also investigate an anisotropic model and show that the 
coupling between planes is essential for the validity of mean-field-type order parameters. 



PACS numbers: 71.10.Fd, 75.10.Lp 

The Hubbard modelS is one of the most important 
and most prominent theoretical models in modern con- 
densed matter physics. Originally introduced in or- 
der to describe magnetism in transition metals, the 
Hubbard model has been most intensively investigated 
since Anderson's proposal that the model should capture 
the essential physics of cuprate superconductors .0 The 
Hubbard-Hamiltonian 
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describes itinerant electrons with spin a on a lattice 
with nearest neighbor sites (ij), interacting through 
short-ranged Coulomb repulsion U. The success of this 
relatively simple lattice model for mobile interacting 
electrons is based on its ability to explain a number 
of important phenomena in condensed matter physics. 
Among these are the (Mott-Hubbard) metal-insulator 
transitionJjQ antifcrromagnetism,El fcrromagnetismja in- 
commensurate phases, u phase separation, Ha and normal- 
state properties of high-T c materials.BE3 

In spite of the prominence of the Hubbard model in 
condensed matter theory and the deceptively simple two- 
parameter form of its Hamiltonian, comparatively little 
is known exactly (or even accurately) about its solution. 

In d = 1 the exact solution has been determined by 
the Bethe Ansatz techniqueO Even in the extreme weak 
coupling regime the ground state of the one-dimensional 
Hubbard model is nonperturbative, possessing a discon- 
tinuity at U = + where the Mott-Hubbard gap opens. 

Much has recently been learned about the solution in 
high spatial dimensions (d — oo), but even here&LL-J the 
exact analytical solution is unknojApL so that one, has, 
to take recourse to perturbativeBEJ~0 or numericaEJH 
methods. There is no doubt that the half-filled Hub- 
bard model exhibits long-ranged antiferromagnetism in 
the high-dimensional limit, whereas the situation away 
from half-filling is less clear. . 

In d = 2 much is known from exact diagcpalizationO or 
Monte Carlo simulation of finite systems.Ej In addition, 
much effort is spent on the analytical3 and numericalcIrE-3 
analysis of perturbative and nonperturbative results. 
While there seems to be now a consensus that at half- 
filling the ground state has long-range antifcrromagnetic 



order, there are still controversies when the system is 
doped and - even at half-filling - concerning the ques- 
tion of whether there is or is not a precursor pseudogap. 

In so far as numerical work on the_<itcee- dimensional 
Hubbard model has been carried out£3t3 the extrapola- 
tion of the results to the thermodynamic limit is clearly 
made difficult by the small linear system size. Much 
effort has been invested in the analysis of the phase 
diagram using perturbative techniques. In the strong 
coupling regime (U — > oo), where the half-filled Hub- 
bard model reduces to an effective Heisenberg model, 
high-temperature series expansions and l/.S'-expansions 
can reliably be used to estimate the Neel temperature 
and the ground state energy, respectively. However the 
analysis of the intermediate and weak coupling regime 
is much more difficult: For this regime a variety pf 
approaches has-been proposed, e.g., by Kakchashi,E3 
Logan,E2l Cyrot,H and Dare.E3 These approaches all yield 
estimates for the Neel temperature which reduce to the 
Hartree-Fock result at weak coupling. However, it is well- 
known from studies of the Hubbard model, baaed on the 
1/d-expansionEj or the local approximation^ that the 
Hartree-Fock approximation overestimates Neel temper- 
ature and order parameter in d — 3 by a factor of the 
order of four, even in the extreme weak-coupling limit. 
The precise value of this renormalization factor in d = 3 
is as yet unknown. 

The goal of this paper is, first, to present exact results 
for the broken-symmetry phase of the three-dimensional 
Hubbard model at weak coupling. In particular, we 
present asymptotically exact formulas for the Neel tem- 
perature and order parameter. It will be seen that this 
asymptotic formula also yields a useful approximation 
formula for the Neel temperature for all U < D, where 
D = 12t is the band width of a three-dimensional hyper- 
cubic lattice. Second, we address sc perturbation theory 
for the broken-symmetry phase in low dimensions (d = 2) 
in the ground state. Our main result here is that sc per- 
turbation theory breaks down altogether, and that low- 
order sc perturbation theory bears no relevance to the 
behavior of the exact solution of the Hubbard model in 
d = 2. We conclude that the antiferromagnetic order pa- 
rameter cannot have a mean-field form. This result has 
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obvious relevance for theories of high-T c superconductiv- 
ity, many of which are based on perturbation theory in 
strictly two-dimensional systems. Third, we investigate 
the order parameter of the anisotropic Hubbard model, 
where two-dimensional planes are weakly coupled in c- 
direction by a small hopping amplitude t±_ -C tn. We 
demonstrate that the sc perturbation series converges - 
however sluggishly - as long as the interplane hopping t± 
is finite. 

Since the behavior of the Hubbard model even in the 
weak-coupling limit is nonperturbative one has to ap- 
ply self-consistent theories. .-There are several ways of 
imposing self-cansistency.EiTLJ The method at fixed or- 
der napimeteiEa and some of its results are described in 
RefsLrEj for the special case of the Hubbard model in 
high spatial dimensions (d — > oo). Here we extend these 
investLgatioas to all finite dimensions. As described in 
RefsJj't3E30 the order parameter A and the Neel tem- 
perature Tn are determined by the roots of the optimiza- 
tion equation 



df „ 
— = 
dA ' 



(2) 



where / is the free energy density. It then follows that 
one has to determine the order parameter from the self- 
consistency condition 
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In Eq. (|[) ho denotes the symmetry-breaking field, 77(e) = 
sgn(e)-\/e 2 + the dispersion, and (3 = l/T the inverse 
temperature. Nd(y) being the density of states in d di- 
mensions. It is important to note that in each order ho 
and A have to be determined self-consistently from the 
Hartree-Fock contribution and the fluctuations together 
to obtain systematic corrections to mean-field theory. 

First we consider sc perturbation theory in the broken- 
symmetry phase in dimensions d > 3. Fortunately, al- 
though the evaluation of the various higher order dia- 
grams in this approach at finite values of the interaction 
(U > 0) is difficult, if possible at all, the results at weak 
coupling (i.e., for U — ► 0) are simple: For all d > 3 one 
finds that the exact value of the Neel temperature Tn 
and the exact order parameter A, can be expressed in 
terms of their Hartree equivalents and a scaling factor 
qd- As is well-known, the Hartree Neel temperature is 
exponentially small for — > 0: 
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where I d can be expressed in terms of an integral, 

Now, we find that the exact expressions for Tn and 
A(T) differ from the mean- field results by a scaling fac- 
tor q d : 
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For the special case d = 00 we know from Ref. [l5| 
that goo = exp(— 7), where 7 = v / 2arcoth(v / 2), so that 
goo ~ 0.2875. The renormalization factor q d for all d > 3 
is determined by the second order diagrams only and can 
be conveniently written in the form of an exponential of 
a lattice sum: 
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where Sd is given by 
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and Fj and Gj(r) are integrals over Bessel functions: 
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This expression is formally exact for all d > 2. We have 
numerically evaluated the renormalization factor qd for 
d — 3 and find the following result: 

q 3 = 0.2599 . 

Thus, the exact Neel temperature and the exact order 
parameter are smaller than the Hartree-Fock predictions 
by almost a factor of four. 

Several remarks are in order: (i) The term labeled 
by j in the lattice sum stems from second order di- 
agrams containing two Hubbard vertices a distance |j| 
apart, (ii) Since_we know from previous calculations in 
high dimensionaij that sites at a relative distance | j | give 
a contribution of order d~^ to the free energy, one ex- 
pects that the lattice sum Sd converges extremely rapidly 
if d is large. (Hi) Keeping only the j = tetm in Sd cor- 
responds to the local approximation of Refo. In d = 3 
we find that g^ oc = 0.2673. Comparison of g^ oc ~ 0.2673 
with the exact result q 3 ~ 0.2599 shows that the local 
approximation works very well even in d = 3. (iv) Note 
that odd lattice sites (i.e., interactions between the two 
sublattices) do not contribute to q d to leading order in 
U for U — > 0. This is due to a different symmetry of the 
Green functions for odd and even values of |j|. (v) Also 
note that our result (||) for the renormalization factor 
breaks down for d = 2 since in this case the integrals Fj 
diverge logarithmically for all even | j | . 
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FIG. 1. Neel temperature Tn as function of U from various 
approaches: Hartree-Fock theory (dot-dashed line), Heisen- 
berg limit from high temperature series expansions (dashed 
line), quantum Monte Carlo simulations in infinite spatial 
dimensionstU (open circles) and three dimensions^ (filled cir- 
cles, the dots are meant as a guide to the eye only), and the 
result of this work (solid line). 

The asymptotic result for the Neel temperature in Eq. 
(Q) has a relevance beyond the pure weak coupling limit. 
This is particularly clear in infinite dimensions , where 
accurate quantum Monte Carlo data are availableJia A 
comparison of the Monte Carlo results with the appro- 
priate asymptotic formula for d = octLa reveals excellent 
agreement for all U < D. Therefore one expects that 
([^) yields an equally good approximation for U < D 
in d = 3. The phase diagram of the three-dimensional 
Hubbard model was calulated with a quantum Monte 
Carlo technique by HirschcJ and Scalettar et AjX^ A 
comparison of the Monte Carlo data from Ref.E3 with 
the expected behavior (Q) reveals a significant discren- 
ancy, suggesting (as was also pointed out in RefsoS) 
that the existing Monte Carlo estimates for Tn are too 
high. Improved Monte Carlo simulations^ show a sig- 
nificantly reduced Neel temperature. Yet it is still more 
than 15% higher than the asymptotic value (see Fig. [I]). 
Our weak-coupling approximation formula (Q) and the 
strong-counling approximation formula = 3.83t 2 /U 
from RefH can hence serve as benchmarks for future 
Monte Carlo simulations. 

Since the upper critical dimension of the Hubbard 
model is presumably d u = 4, one expects that the critical 
behavior of the order parameter in d = 3 is characterized 
by a nontrivial critical exponent for all U > 0. Yet the 
exact asymptotic behavior of A(T) in (Q) displays mean- 
field behavior near Tn, with a critical exponent /3 = 1/2. 
These two observations can easily be reconciled by cal- 
culating the size AT of, the Ginzburg region in d = 3, 
which requires as inputEa the correlation length £ and the 
jump in the specific heat AC at Tn- An explicit calcula- 
tion shows that £ is exponentially large at weak coupling, 

1/2 

£ oc [Tn(Tn — T)\ , while AC is exponentially small, 
AC cx Tn. Combination gives for the relative size of the 
Ginzburg region: AT/Tn oc T n , which is exceedingly 



small at weak-coupling and vanishes for U — * 0. 

Now we address sc perturbation theory for the brpken- 
symmetry phase in d = 2. Since Tn — in d — 2,c2l we 
focus on the renormalization of the ground state order 
parameter: A = q2A H . Calculating only the second or- 
der contribution to Eq. (g) we find d/2/dA ~ hox 2 /6 in 

the limit U -> with x = ^ [in (^J . The solution 

of the self-consistency equation (j^) shows an interaction- 
dependent renormalizationO 
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In contrast to the result in d > 3, where qd is constant, 
the renormalization factor in d = 2 vanishes exponen- 
tially as U — > 0. This result already indicates that the 
sc perturbation series in two dimensions has no small 
expansion parameter at weak coupling. Hence higher 
order fluctuation terms are important to decide whether 
the mean-field prediction for the form of the order pa- 
rameter is at least qualitatively correct or does totally 
fail in d < 2 in the symmetry-broken phase. 

To gain some insight into sc perturbation theory in 
higher orders we calculated bubble and ladder diagrams. 
We find 
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with the function 



F B l(x) =6 + . 



x 
~3~ 



In 



for < x < 1 . The important result is that the optimiza- 
tion condition Eq. (^J) has no roots for < x < 1, apart 
from the high temperature solution A = 0. 

We comment on the results, (i) Our analysis shows 
that Hartree-Fock theory bears no relevance to the be- 
havior of the Hubbard model in d = 2 at half-filling. The 
mean-field result is completely destroyed by quantum 
fluctuations. Hence the antiferromagnetic order param- 
eter must have a completely different form: ln(4t//io) <C 
y/t/U. (ii) A calculation of bubble and ladder dia- 
grams in d = 00 shows that all contributions /„ for 
n > 2 yield small corrections to q of the form q — 
qoc exp [2A?' oo (0)7 2 [/ + 0(U 2 )] , where 7 has been defined 
above, (iii) Concerning the behavior of sc perturba- 
tion theory we conclude that only in dimensions d > 2 
low-order sc perturbation theory is capable of describing 
the physics of the Hubbard model correctly. In low di- 
mensions sc perturbation theory in the symmetry broken 
phase diverges, (iv) Remarkably, in d — 1 the structure 
of sc perturbation theory is exactly the same as in d = 2, 

now with x — ^ In ( ||J in Eq. (^). 

The observation that sc second order perturbation the- 
ory gives even quantitatively exact results in d = 3 
and diverges in d = 2 indicates that between two and 
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three spatial dimensions must be a transition where 
sc perturbation theory becomes inadequate. We intro- 
duced a weak hopping amplitude t± in c-direction. The 
anisotropy is then given by the dimensionless quantity 
A = tx/tu -C 1 where tu denotes the hopping in the 
planes. In this anisotropic three-dimensional model the 
renormalization factor can again be expressed as a lattice 
sum according to Eq. (^|). However, the smaller A is the 
more sites have to be summed over. In the j\ , jVplane all 
points within a region 1/A 2 around the origin contribute 
with ln(l/A) to the renormalization factor, yielding 



qx = exp 



const. 



A 2 ln(l/A), 

The renormalization factor is thus strongly reduced in 
less than three dimensions and zero in d = 2. This result 
demonstrates that the coupling to the third dimension is 
essential for mean-field theories to be qualitatively valid. 

In summary, we have investigated the half-filled Hub- 
bard model in two and three dimensions at weak cou- 
pling calculating systematic corrections to Hartree-Fock 
mean-field theory. We have shown that the mean-field 
antiferromagnetic Hartree-Fock solution does not yield 
exact results in the weak coupling regime. In particu- 
lar, we showed that quantum fluctuation induced correc- 
tions make sc perturbation theory diverge in d < 2. In 
d > 2 the exact order parameter and Neel temperture 
differ from the mean-field result by a scaling factor. 

We acknowledge helpful discussions with F. Gebhard, 
M. Jarrell, A. Muramatsu, and T.M. Rice. 
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